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Abstract 

We compute the energy density at arbitrary temperature of the half 
plane Ising lattice with a boundary magnetic field Hb at a distance M rows 
from the boundary and compare limiting cases of the exact expression with 
recent calculations at T = Tc done by means of discrete complex analysis 
methods. 



1 Introduction 

Recently, Hongler and Smirnov [T] have studied the energy density of the Ising 
model at the critical temperature Tc on a half plane lattice a distance M rows 
from the boundary where two special cases of boundary conditions were consid- 
ered 1) free and 2) fixed with all spins up. The problem was considered on an 
isotropic lattice with a boundary of arbitrary shape. When specialized to the 
half plane the results [2 are for free boundary conditions 

(o-M.oc^M-i.o) - ")buik = -tAt + o{M-^) (1) 

ZttIxI 

and for fixed plus spin boundary conditions 

(cta./.O^TM-1,o) - ">bulk = T^^r + 0[M~^) (2) 

attM 

where in the isotropic lattice the vertical bulk energy density (£")buik is 

('^^^)bulk = (o'Af,oO'A/-l,o)bulk = (3) 

The computations of [li are done by means of discrete complex analysis. 

It is the purpose of this note to compute and study the energy density 
(cM.oCAf-i^o) for the anisotropic lattice at arbitrary temperature on a half plane 
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with a magnetic field Hb applied to the boundary row. The energy operator is 
thus 

oo oo oo 

^ = -^ ^ {Eiaj,k(Tj,k+i + E2(Jj,kCrj+i,k} ~ Hb ^ ai^k (4) 

j — l fc— — oo k— — oo 

where we follow the notations of [3] and let aj^k specify the spin in row j and 
column k. This reduces to the half plane case of [T] when Ei = E2 and T = Tc 
with Hb = for free boundary conditions and Hb — 00 for plus spin boundary 
conditions. The exact result as calculated by Pfafhan methods is given in sec. 
2. Limiting cases with M ^ 00 ioi T < Tc and T Tc— are obtained in sec. 
3 and for T > Tc and T ^ Tc+ in sec. 4. We conclude in sec. 5 with a short 
discussion. 



2 The energy density {(TMflaM-ifi) for arbitrary 
T and Hb 

The computation of (ctm.oO'm-i.o) is done straightforwardly by means of Pfaffian 
methods. The details are given in chapter 7 of [3] and using the result (3.16d) 
on page 152 we immediately find that 



with 



~ 27r 



de- 



(o'M-l,oO'J\/,o) — (^")bulk — I 

1 [{1 - zf) - Z2a{l + zf + 2ziC0s( 



(5) 



a 



2M-1 



[(1 + - 2ai cos6')(l + - 2a2 cos 6*)] 

'2 



1/2 



{e'*^ -l)/(e'<^ + \) + iz^z^^v/v' 



(e^" - l)/(e*« + 1) - iz^Z2^v'/v 



where we use the definitions 



ta.n\\E.j/kBT, 



ai 



Zi(l-Z2) 
(I + 22) ' 



z = taxih Hb/ksT 
a2 



(1-^2) 
Zl(l + Z2)' 



(6) 

(7) 
(8) 



the quantity v/v' is given by (3. 7), (3. 14) and (3.20) on pp 120-122 of [3^, as 



v/v' = 



z|(l + zl + 2zi cos 6*) - Z2(l - zf)a 



(9) 



2ziZ2 sin0 

a is the largest root of (3.2) on page 86 of f3| 

{I + zf){l + zl) - 221(1 - z^)cos6i- Z2(l - z^)(a + a"^) = (10) 
which is explicitly given as 

2.2(1-.^) {a + ^f)i^ + 4)-2ziil-zj)cos9 

+ zi(l + Z2)^ [(1 + al- 2ai cos6')(l + - 2^2 cos^)] ^^^j 

(11) 



2 



where the square root is defined to be positive for real 9. The vertical energy 
density )buik is given by the equivalent forms 

1 r {\-zl)-zMl + zl + 2z,cos6) 

\t ;bulk = Z2 + — I m- 



2ttJ_^ Z2(l-22)(l-a2) 



1 f^^ 



(1 -a3e*^)(l -a^e-'^) 
1 - Q;3e-*^)(1 - a4e*'') 



1/2 



(12) 



with 



^2(1 - zi) (1 - ^1) ,,„x 
"3 — r- , "4 = —77- r (13) 

l + Zl Z2[l+Zl) 

where to obtain the last line of (IT2|) we have used identities of [4]. 
We note in particular that when Hb = 0, 

I=^rd0 i , [(1 f 2^.(1 + ^^2., cos ^)] 

2^y_^ a2A^ 1 [(1 + ^2- 2ai cos 0)(l + a2_ 2^2 cos 0)]'/' 

3 Expansions for M — > oo for T <Tc 

We make contact with the computations of by computing the large M be- 
havior of / as given by © in several cases. 

3.1 T <Tc and Hb = 

When Hb ~ and T < Tc is fixed we obtain the large M behavior of / by 
expanding ([14]) by steepest descents. The maximum of the integrand of (|14p is 
at = and thus expanding for small 6 

lna 6l ~ln — T + 7i ' (^^) 

V (I-21) / (1 - ai)(l - q;2) 

(1 - z2) - z2Q;(1 + zf + 2zi cos 6) 

' ^ {z|(l + zi)2 - (1 - zi)'} + 0(^'), (16) 



1-Zl 

{(1 + q2 _ 2q:iCos6')(1 + - 2 cosa26')}^^^ 

^ {z|(l + zi)2-(l-zi)n + O(02) (17) 



2:1(1 + 22) 
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we find 



_ 1 Z2(l-z2)(l + Zi)2 r 2Af l„a(e) /;Lg^ 

2^ (1-Zi)2 ^'^^ 



Then in (UHl) set 



x^ = 2M^2^ fi^^ ^ (19) 

(1 - ai)(l - 02) 
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and (.\JM — > oo we obtain the result that as M — > cjo 



(o'A/,oO'A/-l,o)^(^^)bulk — — 



Z2(l-z|)(l + zi)2 (l-ai)(l-a2) 



2(1 



2-KZxaiM 



22(1 + 2:1) 



(l-^i) . 
(20) 



-2M 



3.2 T <Tr and > 



The large M behavior of / for T < Tc and Hi, > fixed is obtained from © also 
by steepest descents where in addition to p3|) -([T7 1) we also need the expansion 
for 6* ~ 



(e'» - l)/(e*^ + 1) + iz^z^'v/v' z|A[(l - a2) - z^A] 



where 



A = (1-0^2) 
Thus we find as M 00 



^2(1 + 21)^ 



4z2(l-a2)2 
4 



(1-ai) (l + Z2)2(f-ai) 



z|(l-Z2)(l + 2l)' [(l-«2) 



z2A](l-ai)i/2 



8(1 + Z2)(1 

This is negative for 



zi)^ z2^2^(l-a2)(2ia2A/)3/2 



Z2(l +2l) 



(1-zi) 



< (1 - a2)/A = -(1 + Z2)'(l - ai)(l - ^2) 



and positive for 



z^ > -(l + Z2)'(l-ai)(l-a2) 



(21) 



(22) 



-2M 

(23) 



(24) 



(25) 



We note that both (|20p and have the same exponential decay but that 
([23)) decays faster by a factor of 1/AI than does ([20]) . We also note as z — >■ 
that the amplitude of (1^^ diverges. Therefore in order to connect together the 
regimes of Hh = and Hi, > a crossover regime is required. 

3.3 The crossover regime T < T^, Hf, with z'^M fixed 

The crossover from Hi, — to Hi, > ioi T < Tc is obtained by considering 
z ^ and M ^ 00 with z^M = 0{1). Then when M6^ = 0(1) we have the 



expansion which replaces ([21 



{e'^ -l)/{e'^ + l)+iz^Z2^v/v' 



M6i2 



l)/{e'S + 1) -iz^z^\'/v M6'2 + 4z2M(l-a2)/z|^ 



(26) 
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Then using ((11]) -(H?]),!!!!]) and setting 



(1 - ai)z^A 
we find from ((6]) tliat tlie crossover function is 

{o'Mfi'^M-1,0) — ('S^'')bulk 



(27) 



27r(l-zi)2 



2zia2M 



Z2{l-zl){l + z^f ({l-a^){l-a2)V'^ (Z2{l + z,) 



l-zi 



-2M 



du- 



+ ^2 



When C ^ ^ reduces to (UHl) and when C ^ 00 (QSl) reduces to ^ with 
z ^ 0. 



(28) 



3.4 T ^ Tc- and //f, = 

In order to obtain the resuhs ([1]) and ([2]) of 1 at T = T^, where 

ZlcZ2c + Zlc + Z2c = 1 ^2 = 1 



(29) 



we consider 012 — >■ 1 in the asymptotic expansions for M — >■ cxo (P0)) . (P51) . (|28p 
and observe that the exponential factor 



^2(1 + ^ g4zie(l-a2)/(l-Zle) 

1-zi 



(30) 



and that the coefficients either diverge or vanish. Therefore when 1 — a2 H> a 
separate expansion is needed. In the integral (|6]) for / we set 

0/(1 - 02) = 2; (31) 

with — > and 0:2 —> 1 with x fixed of order one and use the approximations 

2zic 



a ~ 1 



l~zi 



-(l-a2)Vr+^ 



and 



a2Q:[(l — z\) — Z2a{l + zf + 2zi cos6')] 
[(1 + al- 2ai cose')(l + aj - 2 cos6l)]i/2 

Then, defining for M — 00 and 1 — a2 — J' 

4zic 



22:1c 
1-zl 



1 + 



1 - 2^ 



M(l-a2) 



(32) 



(33) 



(34) 



we obtain the result 



(o'A/,oO'A/-l,o) — (^")bulk 
zUl-^2, , 



(35) 
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which, using x = sinhy is rewritten as 

zic(l - a2) 



2zi(l-a2) 



41 -^L) J- 



7r(l-z2) 



[Ki{m) + Koim)] 



dy[coshy + l]e~""'°'^y 

(36) 



where Kn{z) is the modified Bessel function of order n [S]. 
When m — 7> 00 we use the first term in the expansion 



" 2m 



1 + — + 0(m-2) 

8m 



to find that (l36l) reduces to 



^ \ Ip-fX / 2zic(l - 02) ^-4zieAf(l-a2)/(l-g?J 



(37) 



(38) 



which agrees with (|20|) in the limit 02 — J> 1 when we use ([30 
When m — )■ we use 



iiri(m) ~ 1/m, Ko{m) 
to find that (l36l) reduces to 



{<^M,0<^M-l,o) — (^")bu: 



Ik 



Inm 



27rAf 



(39) 



(40) 



which agrees with the result of T for free boundary conditions ([TJ. 
3.5 T ^ Tc- and > 

When T — >■ Tc with i/^ > we need the further approximation that by using 
(EH) in dSl) 



v/v' 



, i-VTT. 



and thus 



(e'^-l)/(e"^ + l) + ^z2z2'^/'"' 
(e'^ - l)/(e^« + 1) - iz^za-^w'/w 

Using (|42]) in dH) with (l3T|)-(l34l) we obtain 

zic(i - Q2) r 'VTT^-i 

nil - 



Vl + 2:2 - 1 



(41) 



(42) 



(43) 
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(o'J\/,00'M-l,o) — (i^^^Obulk — 



which setting x = sinh y gives the result 

= ^']f-^f [K,im)-Ko{m)] (44) 

which differs from (j36]) only in the sign of the term with Ki{m). 
When m 00 we use (l37l) in (l44t to find 



(45) 



(O'M.OO'M-I.o) - (f")bulk = ^ , , (46) 

ZttM 

which agrees with the result of T for fixed spin boundary conditions 

3.6 The crossover regime T T^—, Hb with z'^M fixed 

It is of further interest to determine the crossover between the results (I36p and 
pi]) and the specialization to the crossover between the two results (P) and ©. 
When z -i' and M -> 00 with z^M = 0(1) and when M(l - 02) = 0(1), we 
have the expansion which replaces p2|) 



which agrees with with a2 ^ I- 
When m we use (l39l) to find 



(e'^-l)/(e"^ + l) + zz'z2''^/^' mo;^ [VTT^ - 1] - Cg [VTT^ - 1]' 
(e'« - l)/(e*'' + 1) - iz2z2"iu'/i; "^^^ mx2[Vl + a;2 - 1] + Cca;^ 

where is obtained from (|27l) with Q!2 — > 1 as 
Then using gT]) in ([6]) with (|3T])-(l34l) we obtain, 



j_ ^ic(l-«2) /-"^ d:r ma;2 - [n/T+^ - 1] 
which setting x = sinhy gives the result, 



(49) 



zic(l-a2) /"^ , m(coshy + l)-C, 



2 



n 2w . !- (coshy-l)e— (50) 

- ztc) J -00 m(coshy- 1) + Cc 

When ^ 0, ([55]) is recovered, and when Cc — ^ 00 (|44p is recovered. 
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Finally in order to interpolate between the results ([1]) and ^ we let m — > 
in (l49l) and set mx — t to obtain 



(fM^OfM-l.o) ~ )bulk — — 



1 



2nM 



dt 



(51) 



In Figure [l] the integral in (|5ip is evaluated numerically for various C,'^ . The 
integral vanishes at Q — 0.610058 • • •. 




Figure 1: Numerical plot of the integral in ([FT]) . 



4 Expansions for M ^ oo for T > Tc 

The fundamental result ([5]) holds for T > Tc as well as T < Tc- The analysis of 
the various special limiting cases is parallel to T < Tc where now a2 > I and 
(fT5|) ~ (fT7|) are replaced by 

ln«(^) - In + n T (52) 

Vz2(l + zi)/ (l-ai)(a2-l) 

(1 - zf) - Z2a{l + zl + 2zi cos 9) 

4zi(l-^i) .2 

(l + Z2)2(l + zi)(l-ai)(a2-l) 
{(1 + a? - 2ai cos e){l + - 2q2 cos 6')}^/^ 
= {l-ai){a2-l) + 0{9^) (54) 

4.1 T > Te and Hb = 

Using ([52]) -(HH) in ([6]) we find that M ^ oo, 

(o'A/,Of A/-l,o) — ('^^")bulk 



(53) 
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2a2(l + zi)2(l + Z2) V27r2ia2(l - ai)(a2 - l)M3/2 



^2(1 + 21) 

which is to be compared with the corresponding result pO| for T <Tc- 

4.2 T > and //fe > 
When T > T^, and 6I - 

v/v' 

and ((2T|) is replaced by, 



-2M 

(55) 



(l + Z2)'(l-ai)(«2-l) 
2zi6' 



(e"^-l)/(e''^ + l)+zz2z-lt,/t;^ ^ _ z2(^ ^ ^2)2(^ _ _ ;L)2 

(e''» - l)/(e»« + 1) - iz2z2"iz;'/w 



z|[(a2-l) + z2A]^?2 



(56) 



(57) 



Thus as M — ^ 00 

(o'M,Of Af-l,o) — (^")bulk 
2(1-Zi)222 



Zi(q;2 - 1) 



z2(l + zi)2[(a2 - 1) + z^A] V 2^Afa2(l - ai) 



(l-^i) 



-I -2M 



2:2(1 + 2:1) 



(58) 



with A given by ([22]). The result ([58| is positive for all > q in contrast with 
the corresponding result ([23]) for T < Tc which changes sign at z^ = (1 — 02) /A. 

4.3 The crossover regime T > T^, Hb ^ with z'^M fixed 
In this case we find that 



(e*« - l)/(e»» + 1) - iz'^z:^\'/v 



with defined by ^ and 

= 2Af6'2 

thus we obtain the result 



2lQ!2 



(1 - ai){a2 - 1) 



(59) 



(60) 



(o'A/,00'A/-l,o) ^ (^'')bulk 

(l-zi)2(l-2C2) 



2a2(l + 2i)2(l + Z2)2v/27rzia2(l - ai){a2 - l)M3/2 



(I-21) 



22(1 + 21) 



-2Af 

(61) 



which agrees with when <^ — >■ and agrees with the z — >■ limit of 
when ^ ^ 00. 
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4.4 T Tc+ and Hb = 

Approaching T^. from above, ([5^ is replaced by, 

2zic 



a ~ 1 



1 - z2 



(a2-l)Vr+^ 



where now, 



X = 0/(a2 - 1) 

with 9^0 and a2 — J- 1 with x fixed of order one and p3p is replaced by 



a2a[(l — zf) — Z2a(l + zf + 2zi cos 6*)] 
[(1 + - 2aiC0s6')(l + al - 2cos6l)]i/2 

Defining 



2zic 



1 - 



1 



4z 



Ic 



M(a2 - 1) 



(62) 
(63) 

(64) 

(65) 



and setting x — sinh y we find 



(o'Af,oO'A/-l,o) — ('^^'')bulk 

2zic(a2 - 1) 



zic(q;2 - 1) 



^(1 - ^L) 

-i^o (m) + i^i (to)] 



which is to be compared with (pG]). 
When TO — 7> oo, (l66t reduces to 



dy [-1 + coshy] e-'^''°'^y 



(66) 



(o'Af,00'M-l,o) — ('?'')bulk 
' ' 1 



1 



8 V 2TTZic{a2 - 1) M3/2 



„-42i,Af(a2-l)/(l-zL) 



(67) 



which agrees with ([55)) with a2 — >■ 1. 

When TO — >• 0, (|66l) reduces to the resuh ([T]) 

(o'Af,oO'A/-l,o) — (^^)bulk = — 

4.5 T ^ Tc+ and Hf, > 

When T ^ Tc from above, (PT|) is replaced by 

V Vl + a;2 + 1 



27rM 



(68) 



and dm is replaced by 

(e'» - l)/(e*» + 1) - iz'^z:^^v' /v 



VI + a;2 + 1 



(69) 



(70) 
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Thus setting x = sinh y we obtain the result 

('TM.o^M-i^o) - (f">buik = ^'^l^'^^^ r Mcoshy + l)e-^-oshy 

^^-^^f^^[K,{m) + K,{m)] (71) 

which is to be compared with 
When m — > oo, (|7T|) reduces to 



which agrees with (|58l) with a2 1. 
When m — )■ 0, ([7T|) reduces to 

(ctm.octm-i.o) ^ (^^)buik = " (73) 

zttM 

which agrees with the result ©. 

4.6 The crossover regime T — )• Tc+, Hh with 2;^M fixed 
In this case we have 
(e'f - l)/(e*^ + 1) + tz^z^\/v' ffix^iVTT^ + I] - CciVTT^ + 1]^ 



(e'^ - l)/(e*'' + 1) - iz^z^^v'/v mx^[Vl + x^ + I] + Qx"^ 

Thus, using (j64l) and setting a; — sinhy, we obtain 

(o'j\/,oO'A/-l,o) — (^")bulk 

zic{a2 - 1) /""^ , m(coshy - 1) - Cc 



(74) 



7r(l-zL) 



which is to be compared with (|50p . When — ^ 0, ([66|) is recovered, and when 
— > oo ([71]) is recovered. When rri — )■ the result (fSTj) is again obtained. 

5 Discussion 

In this paper we have derived leading behavior of the energy density operator of 
the Ising model on an anisotropic lattice M rows from a half plane boundary at 
critically with a magnetic field Hb on the boundary by use of Pfaffian methods. 
When the field is zero and infinity we regain the results ([T]) and ^ obtained in 
[1] by means of discrete complex analysis. Furthermore in (jSip we have obtained 
the result in the more general situation where H^M is fixed with Hi, ~> and 
M — >■ oo. This result goes beyond the computations of [T| and we have obtained 
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many results for T ^ Tc- It is of interest to obtain these results also by the 
methods of discrete complex analysis. 

We would also like to take this opportunity to remark that it would be most 
useful to extend the results of [T] in several directions. One such direction is to 
consider discrete complex analysis on surfaces of higher genus and to derive and 
extend the results of and 

A second direction is to consider inhomogeneous random lattices. One such 
case is the layered Ising model where the vertical interaction constants are the 
same in all columns but are chosen randomly from row to row [31 chapters 14 
and 15] where it is known that there is an entire temperature region around 
Tc where the correlation functions are algebraic. All of these problems can be 
considered as problems with free fermions and thus the methods of discrete 
complex analysis should apply. 
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